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Abstract 

The effect of wall cooling on the two-dimensional linear stability of subsonic flows 
over two-dimensional surface imperfections is investigated. Results are presented for 
flows over smooth humps and backward-facing steps with Mach numbers up to 0.8. 
The results show that, whereas cooling decreases the viscous instability, it increases 
the shear-layer instability and hence it increases the growth rates in the separation 
region. The coexistence of more than one instability mechanism makes a certain 
degree of wall cooling most effective. For the Mach numbers 0.5 and 0.8, the 
optimum wall temperatures are about 80% and 60% of the adiabatic wall 
temperature, respectively. Increasing the Mach number decreases the effectiveness 
of cooling slightly and reduces the optimum wall temperature. 


I. Introduction 


Due to proven achievability of Natural Laminar Flow' (NLF), there is an increasing 
interest to use it for the design of high performance aircraft. The substantial drag 
reduction with NLF has promoted more analyses of ways to achieve and maintain NLF 
on airfoils and other aerodynamic geometries. The maintenance of NLF is critically 
sensitive to the location of transition, which is strongly affected by surface 


imperfections. Since many of these imperfections cannot be avoided by modern 
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manufacturing techniques 2 3 , a guide is needed for their allowable sizes and methods 
that should be used to control their effect on NLF. 

There are some empirically based criteria in the literature 45 for prediction of the 
transition of flows around imperfections. But these criteria are for special cases and 
geometries, and they do not explain the instability mechanisms enhanced by the 
imperfections or the physics of ways to control them. Nayfeh, Ragab, and 
Al-Maaitah' studied analytically the stability of incompressible flows around 
two-dimensional (2-D) bulges; their method accounted for both viscous and 
shear-layer instabilities. They also correlated their results with the experiments of 
Walker and Greening 7 and found that the e N method can be used to predict transition 
around 2-D bulges. Bestek et al 8 solved the unsteady Navier-Stokes equations using 
finite-differences; they concluded that the unsteadiness of the separated flow can be 
regarded as a phenomenon governed by the hydrodynamic instability. Their 
calculations show that Tollmien-Schlichting waves amplify considerably once the 
separation is enhanced in agreement with the results of Nayfeh et a I* . Burnel et al 9 
aid Gougat and Martin' 0 experimentally investigated the flow over 2-D imperfections. 
Their measurements of the amplified disturbances show that they damp down in the 
region of favorable pressure gradient. A similar trend is seen in the calculations of 
Nayfeh et al*. 

Recently, the effect of compressibility on the achievability of NLF has received 
more attention for non-lifting surfaces. High subsonic and supersonic Mach-number 
flows develop density gradients across the boundary layer, which provide additional 
damping to 2-D and axisymmetric T-S waves. For certain geometricies this 
advantage can be offset by the increase in the adverse pressure gradients. Vijgen 
et al” showed that increasing the freestream Mach number has a stabilizing effect 
on subsonic laminar boundary layers over fuselages. Their Mach number varied from 
low subsonic to 0.8. Hastings et al' 2 reported that NLF extended as far as 37% on a 
NLF fairing installed on a turbo fan nozzle. 

In spite of the previous investigations more understanding of the physics of the 
instability of such flows and ways to control them are still needed, especially for 
compressible flows. The two most common ways for laminar flow control are wall 
cooling in air (or heating in water) and wall suction. Wall cooling stabilizes 
incompressible flows over flat plates in air and destabilizes them in water' 3 -' 0 . 
Mack' 7 and Malik' 8 found that for compressible flows, wall cooling stabilizes the first 
mode but destabilizes the second mode. The question to be asked is how does wall 


V 


2 


cooling affect flows around surface imperfections, and whether these flows can be 
stabilized by this technique, especially, since there is an interaction of more than one 
instability mechanism. 

The purpose of this work is to study the effect of wall cooling on the subsonic 
two-dimensional stability of boundary layers around two-dimensional smooth 
backward-facing steps and humps. 


II. Mean Flow 


The sizes of the two-dimensional imperfections under consideration are such that 
strong viscous-inviscid interactions and small separation bubbles are unavoidable. 
The conventional laminar boundary-layer formulation cannot predict such flows. An 
alternative is to solve the full Navier-Stokes equations, but in such a case the grid 
should be >o fine that important flow characteristics are not smeared by the 
truncation error and artificial dissipation. However, due to the large number of cases 
need to be investigated, solving the full Navier-Stokes equations is a very expensive 
task. A more economical alternative is to solve the interacting boundary-layer (IBL) 
equations or the nonlinear triple-deck equations. 

We calculated the two-dimensional compressible laminar boundary layers over 
flat plates with 2-D surface imperfections using the interacting boundary-layer 
equations (IBL). The flowfield is assumed to be governed by the steady compressible 
boundary-layer equations as follows: 


x-momentum equation 


du , du 
pU ~dx +pV ~dy 


dp i d_( du \ 
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continuity equation 
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and equation of state for perfect gas 

pT = P e T e 


( 3 ) 


(4) 


where velocities are normalized with respect to the freestream velocity IT,, lengths 
are normalized with respect to L*, which is the distance from the leading edge to the 
center of the imperfection, and the temperature and the viscosity and 
thermal-conductivity coefficients are normalized with respect to their freestream 
values T^, , and k„, respectively. Here, 


Re = 



Pr = 



K 


on 



(5) 


where C‘ and Ct are the gas specific-heat coefficients at constant pressure and 
volume, respectively. 


The boundary conditions at the wall are 

u = v = 0 and T = T w at y = f[C(x)] (6) 

where f(Q is the shape of the wall with the imperfections. For the step 

KO = \ htf + erf(0], C = Re~ 3ls X 5l \x - 1) (7) 


and for the hump 


f( 0 = 


h{ 1 - 12C 2 + 16C 3 ), if c < 0.5 
0.0, if ( > 0.5 


where x w is the width of the hump, erf is the error function, h is the height of the step 
or the hump, and X = 0.332057. Away from the wall 

u — > U e and T -* T e as y -+ oo (8) 

where the subscript e stands for edge variables. 
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Using the Prandtl transposition theorem, 


z = y- /[£(*)], w = v - u- 


we rewrite Eqs. (1)-(3), (6), and (8) as 

du , du 



dp 1 

8 

dz 

dx Re 

dz 

d(pw) 

, d(pw) 

0 




u = w = 0,T=T w at z — 0 
U —* U e and T -> T e as z -* oo 

Next we use the Levy-Lees variables 


«*) = f ; 


P e H e U e dx and ti(x,z) = 


jReUe r 

Jo' 


and transform Eqs. (10)-(12) into 


2 t FF ! + VF >i-if( e -jif) + Po(F 2 -Q) = 0 

2 Iff* + V n + F — 0 

2Wt + w ,-J-(^i2.) 

2 U l 2 

-(y- i)/wi-^-< = o 


F = O = -I- 
"e ’ T e 

[21 

v ~ 77J7 7 ['^ r '” v + 
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(19c) 


0 = and /? 0 = 2< ° ^ 


PeUe 


U e dt 


The boundary conditions become 

F — V - 0 and Q = Q w at rj — 0 
F — ♦ 1 and Q — ► 1 as tj —* oo 

F = F(to, *l) and Q = Q(£o. V) at £ = £ 0 


(20a) 

(20b) 

(20c) 


To account for viscous-inviscid interactions, we need to calculate the inviscid 
flow over the displaced surface. This is done through the interaction law, which 
relates the edge velocity to the displacement thickness. Using thin airfoil theory, we 
obtain 




J LE 


d{\np e )/dt 

x-t 


dt + 


1 _ f°° 

Pn J LE 


d{U 9 d)ldt 

x-t 


dt 


( 21 ) 


where /? = Jt — Ml, , the displacement thickness 3 is given by 


i fit 

3 “ JFF J 0 <0 “ W (22) 

and U, is the inviscid surface velocity in the absence of the boundary layer, which, in 
the case of small imperfections, can be expressed as 


u *’ 1+ irl 1 


dfjdt 


pn J LE x-t 


dt 


(23) 


Defining x - f + U e S, we rewrite Eq. (21) as 


U e = 1 + 


pn 


Poo 

L 


u,s 


d{\np e )ldt 

x-t 


dt + 


1 _ f°° 

Pn J LE 


dxjcH 

x — t 


dt 


The principal values of the integrals in Eq. (24) are assumed. 


(24) 


Following Davis and Werle 20 and Nayfeh et al e , we integrate the interaction law 
by parts to eliminate the derivative of x ■ We assume x to vary linearly over a 
differencing interval to obtain a second-order quadratic expression for the edge 
velocity. Furthermore, we calculate the second term in Eq. (24) explicitly from the 
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previous iteration. By assuming the flow far away from the imperfection to be a 
Blasius flow, we can write the interaction law as 

V N + 4>P oi = $ (25) 

For a definition of </> and ^ and a detailed derivation of Eq. (25), we refer the reader 
to Ref. 6. 

Equations (16)-(18) and (20) are solved simultaneously with Eq. (25) using central 
differences in the vertical direction and three-point backward differencing in the 
streamwise direction. 

Ragab, Nayfeh, and Krishna*’ compared the IBL calculations with solutions of the 
thin-layer compressible Navier-Stokes equations obtained using the computer code 
"ARC2D" developed at NASA Ames (Version 1.5 GAMMA). The results obtained using 
the IBL agree very well with those obtained using the Navier-Stokes solver; this is 
true for both the mean flow and the stability characteristics. In Ref. 21 insulated wall 
conditions were used. In the present work, we investigate the influence of continuous 
and strip cooling. 


III. Stability Calculations 

In this work, we consider the linear two-dimensional quasi-parallel stability of the 
mean flow calculated using the interacting boundary-layer formulation. The 
quasi-parallel assumption was justified a posteriori by Nayfeh et al 6 . They found that 
the wavelengths of the disturbances are the order of the boundary-layer thickness. 

The calculations are performed for constant specific heats and Prandtl number. Since 
we are limiting our calculations to subsonic flows, this assumption has a small effect 
on the accuracy of the stability results. Moreover, the viscosity and 
thermal-conductivity coefficients /xand k are assumed to be functions of temperature 
only. Since Pr and C' are constant we take k = h . 

To derive the stability equations, we superimpose 2-D disturbances on the mean 
flow calculated using the interacting boundary-layer formulation to obtain the total 
flow quantities 


p = Pm(y) + p(*.y. 0 

(26a) 

u = u m (y) + u(x,y,t) 

(266) 
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II 

(26c) 

p = Pm(y) + p(x,y,t) 

(26d) 

m = p m (y) + p(X'y.t) 

(26e) 

X = XJy) + X(x,y,t) 

(26f) 

f = T m (y) + T(x,y,t) 

(26g) 


where A and n appear in the definition of the bulk-viscosity coefficient k as 

k = A + - jfi 

The subscript m refers to mean-flow quantities and the overbar refers to total flow 
quantities. Since and A are functions of temperature only, we have 

cfA d 

X = -rf-T=X' m (T m )T and n = T = (27) 

u 1 m u 1 m 

Substituting Eqs. (26) and (27) into the 2-D compressible Navier-Stokes equations, 
subtracting the mean-flow quantities, and linearizing the resulting equations, we 
obtain 
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( 31 ) 


Pm\ 


dT 

dt 


+ u. 
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■ n 1” 


RPr dy dy 


RPr 


dy* 


where 



The linearized equation of state for a perfect gas is 

V M loP = PmT+ pT m 

or 

P = (V^lcP ~ PrrJMTm (35) 

The boundary conditions are 

u = v = 0, T=0 at y = 0 (36) 

u f v,p,T->0 as y -+ oo (37) 

Since the coefficients in Eqs. (28)-(31), (36), and (37) are functions of y only, we 
seek normal-mode solutions of the form 


q — q(y ) exp{/J adx — iiot} + complex conjugate (38) 

where q stands for (u, v, p, T), a is the wave number, and co is the frequency. For 
spatial stability analysis a is complex and c o is real, whereas for temporal-stability 
analysis co is complex and a is real. In this work, we analyze the spatial stability case 
and determine (o from the nondimensional frequency F as co = F/R. 
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Dropping the hat from q for convenience and defining 


Q = (o — au n 


( 40 ) 


we find that a is governed by the eigenvalue problem 


„ DT m iLln j qt 
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and p m = 
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yMl 


(47) 
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IV. Numerical Results 


Following the procedure described in section II, we calculated the mean flow 
over a backward-facing step for different wall temperatures. Figure 1 shows the 
influence of cooling on the skin-friction coefficient 


_ 2Uln(T w ) dF 
f dr > * 


(48) 


In the case of cooling, the separation point is almost fixed but the reattachment point 
moves slightly upstream, resulting in a slightly smaller separation bubble. Moreover, 
cooling causes a larger negative shear prior to reattachment. The pressure 
coefficient C„ = 2{p t T t - 1)/(yA4J„) is plotted in Fig. 2. Although cooling does not have 
much effect on C p far away from the imperfection, it causes steeper adverse and 
favorable pressure gradients around the separation bubble. In Fig. 3, we compare 
the mean profiles of the flows over adiabatic and cooled walls at several locations. 
In general, cooling results in fuller velocity profiles compared to the adiabatic case. 
In the separation region, although the mean-velocity profiles are still fuller away from 
the wall for the cooled wall case, they develop inflection points close to the wall, and 
more negative flows develop near reattachment. The corresponding temperature 
profiles are shown in Fig. 4. The combined effect of cooling on the velocity and 
temperature profiles is a movement of the generalized inflection point closer to the 
wall in the separation region, as shown in Fig. 5. 


For a given mean flow, o>, and R, we solved for the eigenvalue a and the 
eigenfunctions, and then determined the amplification factor from 


N = -\ R 2« t dR (49) 

J *o 

where ft 0 corresponds to Branch I of the neutral stability curve and a, is the imaginary 
part of a The eigenvalue problem was solved using the second-order 
finite-difference subroutine DBVPFD 22 , which is much faster than SUPORT 23 ; the 
results of DBVPFD are in full agreement with those of SUPORT. In all the cases, the 
results are for the most dangerous frequency, defined to be the one that results in 
an N factor of 9.0 in the shortest distance 62 ’. In the case of cooling it is the one 
resulting in the largest N factor in the shortest distance. It turns out that the most 
dangerous frequency for cooled and adiabatic wall are about the same. Figure 6 
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shows the N factor for various frequencies when T m = 0.55 T ad . It appears that 
F = 50x10' produces the largest N factor. 

We analyzed the stability of the mean profiles calculated using the IBL code for 
flows over a backward-facing step. The present analysis accounts for both viscous 
and shear-layer instabilities in the separation region. The effect of wall cooling on the 
stability of such flows is different from that on the stability of flows over flat plates. 
Figure 7 shows the growth rates for the cases of adiabatic and cooled walls. The 
temperature of the cooled wall is 55% of the adiabatic wall temperature T, d . We note 
that the instability is due to the viscous mechanism in the attached region and due 
to a combination of the viscous and shear-layer mechanisms in the separation region. 
Figure 7 shows that cooling decreases the growth rates and hence it is stabilizing in 
the attached flow regions because cooling produces fuller velocity profiles. On the 
other hand, in the separation region cooling increases the growth rates due to the 
increase in the negative shear flow in the separation bubble and the movement of the 
generalized inflection points closer to the cooled wall. However, the growth-rate 
curve corresponding to the cooled case is narrower around the peak value than that 
corresponding to the adiabatic case because the cooled flow reattaches ahead of the 
adiabatic flow. Figure 8 shows the variation of the growth rate with streamwise 
distance for different wall temperatures. Decreasing the wall temperature 
destabilizes the flow in the separation bubble and stabilizes it in the attached flow 
region. Moreover the growth-rate curve gets narrower as the wall temperature 
decreases. Consequently, the overall effect of cooling as measured by the N factor 
depends on the wall temperature, as shown in Fig. 9. When T w = 0.95T id , the 
boundary layer is completely stable ahead of separation and the overall N factor is 
less than that of the adiabatic wall. As T w decreases below 0.8 T, d the increase in the 
growth rates in the separation region overcomes the reduction in the growth rates 
elsewhere, and the net result is an increase in the maximum N factor. 

For a Mach number of 0.8, Fig. 10 shows the growth rates for various wall 
temperatures. The increase in the growth rates in the separation region when M M 
= 0.8 is less than that when M M = 0.5. Moreover, the peak growth rate when 
= 0.8 is wider than that when - 0.5. This makes the optimum wall temperature 
to be at T w = 0.6T ad as it is clear from the resulting N factors shown in Fig. 11. Figure 
1 1 also shows that at M = 0.8 cooling has a slightly smaller effect than at = 0.5 . 

The previous results seem to be general and apply to other imperfections. For 
example, Figure 12 shows the growth rates for a flow at M^ = 0.8 around a cubic 
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hump. The hump width is 0.2 L' and height is 0.003 L*. The behavior of the growth 
rates with cooling is similar to that of the step case. The resulting N factors are 
plotted in Fig. 13. This figure shows that the optimum wall temperature is about 
T w = 0.8 T^. When T w = 0.6T ad the N factor is about the same as that for the adiabatic 
wall. 

In an attempt to lower the growth rates in the separation region, we performed 
calculations for a wall that is cooled everywhere except in the separation bubble. The 
results indicate that this distribution has a very small effect on the N factor as shown 
in Fig. 14. The growth rates are slightly changed as shown in Fig. 15. The reason for 
this small effect is that the mean profiles in the separated region are influenced by 
the cooling ahead of the separation bubble due to the nonsimilarity of the boundary 
layer. 

Next we show the influence of the step position on the N factors for a constant 
frequency. The results shown in Fig. 16 are for = 0.5 and adiabatic wall 
conditions. Unlike the incompressible case the most dangerous step location is not 
the one corresponding to Branch I of the neutral stability curve, but it is the one 
corresponding to a distance half-way between Branches I and II. Figure 17 shows 
that a similar trend is true for the cubic hump. 

From the previous results it is clear that the coexistence of viscous- and 
shear-layer instability mechanisms complicates the effect of cooling on the stability 
of such flows. Since cooling decreases the viscous instability and increases the 
shear-layer instability, there exist an optimum wall temperature that considerably 
reduces the amplification factor. 
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Figure Captions 


Figure 1. Effect of wall cooling on the shear coefficient for a flow over a 
backward-facing step when the step height = 0.003, step slope = - 4.34695, 

M„ = 0.5, Re = 1.0x10', and Pr = 0.72: TJT. d = 1.0,-- TJT, d = 

0.8, and - - - TJT, d = 0.55 . 

Figure 2. Effect of wall cooling on the pressure coefficients for a flow over a 
backward-facing step when the step height = 0.003, step slope = -4.34695, 

A4 00 = 0.5, Re = 1.0x10', and Pr = 0.72: TJT, d = 1.0, TJT, d = 0.8, 

and - - - TJT, d = 0.55 . 

Figure 3. Effect of wa'I cooling ot; the streamwise velocity profiles along the plate 
when the step height = 0.003, step slope = -4.34965, = 0.5, 

Re = 1.0x10*, and Pr = 0.72. The profiles correspond to the following values 
of R starting from left to right: 985, 992, 997, 1012, 1027, 1037, 1042, and 
1051. The separation bubble for the cooled wall starts at R = 997 and ends 
at R = 1037: TJT ad = 1.0, and - - - TJT >d = 0.55. 

Figure 4. Effect of wall cooling on the temperature profiles along the plate when the 
step height = 0.003, step slope = -4.34965, M ao = 0.5, Re = 1.0x10', and Pr 
= 0.72. The profiles correspond to the following values of R starting from 
left to right: 985, 992, 997, 1012, 1027, 1037, 1042, and 1051. The separation 

bubble for the cooled wall starts at R = 997 and ends at R = 1037: 

TJT, d = 1.0, and - -- TJT 3rf = 0.55. 

Figure 5. Variation of the inflection point across the separation region for adiabatic 
and cooled-wall conditions when the step height = 0.003, 

step slope = — 4.34695 , M ^ = 0.5, Re = 1.0x10', and Pr = 0.72: 

T w /T Bd = 0.8, and. — TJT, d = 1.0. 
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Figure 6. The influence of the frequency on the N factor when T w = 0.55T, d , step 
height = 0.003. step slope = - 4.34965, M oe = 0.5, Re — 10x10'* , and Pr = 
0.72: F- 40x10-6, F = 60x10'*, and — .. F = 50x10 *. 


Figure 7. The growth rates for the flow over cooled and adiabatic walls when the step 
height = 0.003, step slope = - 4.34965, = 0.5, Re — 10x10*, Pr = 0.72, 

and F = 50x10 *: TJ T >d = 1 .0 and TJT, d = 0.55. 

Figure 8. Variation of the growth rates with wall temperature: step height = 0.003, 
step slope = - 4.34695, A4„ = 0.5, Re = 1.0x10*, Pr = 0.72, and F = 
50x10-*. 


Figure 9. Variation of the amplification factor with wall temperature: step height = 
0.003, step slope = - 4.34695, M„ — 0.5, Re = 1.0x10*, Pr = 0.72, and F 
= 50x10-*. 


Figure 10. Variation of the growth rates with wall temperature: step height = 0.003, 
step slope = - 4.34695, M „ = 0.8, Re = 1.0x10*, Pr = 0.72, and F = 
50x10-*. 


Figure 11. Variation of the amplification factor with wall temperature: step height = 
0.003, step slope = - 4.34695, = 0.8, Re = 1.0x10', Pr = 0.72, and F 

- 50x10 *. 


Figure 12. Variation of the growth rates with wall temperature for a cubic hump: hump 
height = 0.003, x w = 0.2, = 0.8, Re = 1.0x10*, Pr = 0.72, and F = 50x10 *. 


Figure 13. Variation of the amplification factor with wall temperature for a cubic hump: 
hump height = 0.003, x w = 0.2, = 0.8, Re = 1.0x10*, Pr = 0.72, and F 

= 50x10-*. 
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Figure 14. Effect of an adiabatic separation region on the amplification factor: step 
height = 0.003, step slope = - 4.34695, = 0.5, Re = 1.0x10*, Pr = 0.72, 

and F = 50x10 *: adiabatic wall, — •• TJT t<l = 0.8, - - - cooled 

everywhere except in separation region. 


Figure 15. Effect of an adiabatic separation region on the growth rates: step height 
= 0.003, step slope = - 4.34695, = 0.5, Re = 1.0x10*, Pr = 0.72, and 

F = 50x10 • : adiabatic wall, — •• TJT, d = 0.8, - - - cooled everywhere 

except in separation region. 


Figure 16. Effect of step location on the amplification factor: step height = 0.003, step 
slope = - 4.34695, = 0.5, Re = 1.0x10*. Pr = 0.72, and F = 50x10 *. 


Figure 17. Effect of hump location on the a 
hump width = 0.2, M - 0.5, Re 


mplification fector: bump height = 0.003, 
- 1.0x10*, Pr = 0.72, and F = 50x10'* . 
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